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ABSTRACT 

The prevalence of radio-loud AGN activity in present-day massive halos is determined 
using a sample of 625 nearby groups and clusters selected from the Sloan Digital Sky 
Survey. Brightest group and cluster galaxies (BCGs) are more likely to host a radio- 
loud AGN than other galaxies of the same stellar mass (by below a factor of two 
at a stellar mass of ~ 5 x lO n M0, but rising to over an order of magnitude below 
IO^Mq). The distribution of radio luminosities for BCGs does not depend on mass, 
however, and is similar to that of field galaxies of the same stellar mass. Neither the 
radio-loud fraction nor the radio luminosity distribution of BCGs depends strongly on 
the velocity dispersion of the host cluster. The radio- AGN fraction is also studied as 
a function of distance from the cluster centre. Only within 0.2 r2oo do cluster galaxies 
exhibit an enhanced likelihood of radio-loud AGN activity, which approaches that 
of the BCGs. In contrast to the radio properties, the fraction of galaxies with optical 
emission-line AGN activity is suppressed within r2oo i n groups and clusters, decreasing 
monotonically towards the cluster centre. 

It is argued that the radio-loud AGN properties of both BCGs and non-BCGs can 
naturally be explained if this activity is fuelled by cooling from hot gas surrounding 
the galaxy. Using observational estimates of the mechanical output of the radio jets, 
the time-averaged energy output associated with recurrent radio source activity is 
estimated for all group / cluster galaxies. Within the cooling radius of the cluster, 
the radio-mode heating associated with the BCG dominates over that of all other 
galaxies combined. The scaling between total radio- AGN energy output and cluster 
velocity dispersion is observed to be considerably shallower than the ~ a* scaling 
of the radiative cooling rate. Thus, unless either the mechanical-to-radio luminosity 
ratio or the efficiency of converting AGN mechanical energy into heating increases by 
2-3 orders of magnitude between groups and rich clusters, radio-mode heating will 
not balance radiative cooling in systems of all masses. In groups, radio-AGN heating 
probably overcompensates the radiative cooling losses, and this may account for the 
observed entropy floor in these systems. In the most massive clusters, an additional 
heating process (most likely thermal conduction) may be required to supplement the 
AGN heating. 

Key words: galaxies: active — radio continuum: galaxies — galaxies: clusters: gen- 
eral — cooling flows — X-rays: galaxies: clusters 



1 INTRODUCTION 

Brightest cluster galaxies (hereafter BCGs) occupy a unique 
position at the centre of the gravitational potential well of 
clusters of galaxies. They are amongst the most luminous 
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galaxies in the Universe, up to ten times brighter than typi- 
cal ellip tical galaxies, wit h large characteristic radii (tens of 
kpc; e.g. ISchomberti 1 986) and, in some cases (the 'cD' galax- 
ies; e.g. bemlerlll97^V ~diffusc envelopes which can extend 
hundreds of kpc. The properties of BCGs are found to be 
correlated with those of their surrounding cl usters fe.g. lEdgel 
ll99ll : lLin fc Mohrll2004l ; lBrough et al.ll2005l ), and the forma- 
tion history of BCG and cD galaxies is believed to be distinct 
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from that o f other elliptical galaxies, due to their special lo- 
catio n (e.g. iBernstein fc~ Bhavsar 2001; |p"e Lucia fe Blaizot] 
120061 . and references therein). They are also offset from 
the standard scaling relat ions of early-type galaxies (e.g. 
Ivon der Linden et alj|2006h . 

Brightest cluster galaxies have long been recognised to 
show different radio properties to other cluster galaxies, 
being much more likely to be radio-lu minous than other 
non-dominant cluster ellipticals (e.g. Burns et al.l Il98ll : 
IValentiin fc Biileveldl 11983 : iBurnsi Il990l ). This radio-loud 
AGN activity in BCGs has been proposed as a potential so- 
lution to the cooling-flow problem. Gas in the central regions 
of clusters of galaxies often has radiative cooling timescales 
very much shorter than the Hubble time and, in the ab- 
sence of a heating source, a cooling flow would be expected 
to develop, whereby the temperature in the central regions 
of the cluster drops and gas flows inwards at rates of up to 
~ 1000M Q yr -1 (see lFabianl ll994. for a review). However, 
recent XMM-Newton and Chandra observations of cooling 
flow clusters have shown that the temperature of cluster 
cores does not fall below ~ 30% of that at large radii, and 
that the amount of cooling gas is only a bout 10% of that pre- 
dicted for a classical cooling flow (e.g 
iDavid et alll200ll ; iTamura et ail 2001 



Peterson et al 1120011 : 
Kaastra et al.ll200ll ). 



This implies that some heating source must balance the ra- 
diative cooling losses, preventing the gas from cooling fur- 
ther. 

Heating by radio sources associated with the BCGs 
has gained popularity in recent years, as X-ray observa- 
tions have revealed bubbles and cavities in the hot in- 
tracluster medium of some cluster s, coincident wi t h the 
lobes of the radio sources (e.g. iBohringer et all 11993": 
Carilh et al] Il994l: jMcNamara et al.l l200d : iFabian et al.l 



200d : iBlanton et alJ l200ll ). These are regions where rel- 



ativistic radio plasma has displaced the intracluster gas, 
creating a low-density bubble of material in approximate 
pressure balance with the surrounding medium, which then 
rises buoyantly and expands. Hydrodynamic simulations 
have been able to reprod uce bubbles with properties similar 
to those observed (e.g. Churazov ct a . 20Ql|; iQuilis et al.l 



l200ll ; iBriiggen fc Kaiserll2002l ; iBriiggenllioolh and have also 
shown that, provided that radio source activity is recur- 
rent, the total energy provided by AGN activity can be 
sufficient to balance the cooling radiation losses through re- 
peated prod uction of jets, buoyant bubbles and associated 



shocks (e.g. iDalla Vecchia et al.ll2004l ; IBriiggen et al 



Nusser et al.ll2006l). O ther a uthors (e.g. lOmma et al 



2005 



2004 



Brighenti fc Mathews! 120061 ') have argued that momentum- 



driven jets are an alternative means of distributing AGN 
energy throughout the intracluster medium. 

Very deep Chandra observation s of the Perse u s and 
Virgo clusters (IFabian et al.l 12003" : IFabian et all 120061 : 
iForman et"ai1l2005l . l200dT have revealed the presence of ap- 
proximately spherical weak shock waves in t hese clusters, 
extend ing out to several tens of kpc radius. IFabian et al.l 
l|2003t ) argued that these 'ripples' are excited by the expand- 
ing radio bubbles, and that dissipation of their energy can 
prov ide a quasi-continuous h eating of the X- r ay gas (see 
also iRuszkowski et all 120041 ). iFuiita fc Suzukil i|2005l ) and 
iMathews et al.l "i 2006h argued, however, that if all of the 
wave dissipation occurs at the shock front then too much 
of the heat is deposited in the inner regions of the cluster, 



and not enough out at the cooling radius, le ading to tem- 
eratu re profiles at variance with observations. iNusser et al.l 
2006) suggested that a natural solution to this problem 
is to consider the heating effects of AGN act ivity from all 
cluste r galaxies, not only BCGs. More recently. IFabian et al.l 
(2006) showed that the shocks occurring in the Perseus clus- 
ter are isothermal, meaning that thermal conduction must 
be important in mediating the shock. They argue that this 
prevents the accumulation of hot shocked gas in the inner 
regions, and that the energy of the waves can be dissipated 
at larger radii by viscous effects. 

One critical question for all of these models is whether 
the heating generated by an episode of radio-loud AGN ac- 
tivity is sufficient, during its lifetime, to offset the radia- 
tive cooling losses of the cluster. Viewed in terms of the 
radio bubbles (which are widely considered to be the driv- 
ing force behind the shock and sound waves) , for some clus- 
ters the (pV) energy contained within the evacuated bub- 
bles has been shown to be sufficient to balance the cooling 
losses, at least for a shor t period of time (a f ew x 10 7 yr; 
e.g. IFabian et ail 12003": IBirzan et all |2004 iDunn et all 
120051) . However. iBirzan et all <|2004F showed that this is not 
true for about half of the cluster s they studied (see also 
iRaffertv et aill2006l ). INusser et all l|2006l ) and others have 
argued that the mechanical energy injected into the clus- 
ter by weak shocks may be up to an order of magnitude 
larger than pV, which would ease this problem. An alterna- 
tive possibility is that the heat balance occurs in a quasi- 
static manner, with recurrent low luminosity radio activity 
being punctuated by occasional major eru ptions supplying 
much more energy |Kaiser fc Binnevll20oj ). 

A second important question is what the duty cycle of 
this AGN activity is. The duty cycle determines the rate of 
production of radio bubbles, or equivalently, the timescale 
between the sound wave 'ripples', and hence is required to 
calculate the time-average d heating rate associated with 
AGN activity. iBuriisI (|l990h showed that as many as 70% of 
cD galaxies at the centre of cooling flow clusters are radio- 
loud, but this result was based on a sample of only 14 such 
systems, and not all BCGs are cD galaxies. 

Recently, iBest et all (|2005l ) used data from the 
Sloan Digital Sky Survey (SDSS; lYork et all l2000l : 
IStoughton et al.|[2002l ) to investigate the origin of radio-loud 
AGN activity, and found that the dominant factor which de- 
termined whether or not a galaxy is radio-loud is its mass: 
the fraction of galaxies which were found to host radio- 
loud AGN scaled as M, 2 ' 5 or M^, where A'/* and M B h are 
the stellar and black hole masses of the galaxy, respectively. 
The distribution of radio luminosities was found to be the 
same regardless of the galaxy mass. Combining these re- 
sults with a conversion between radio luminosity and the 
mechanical ene rgy supplied by a j et to its surr o undin gs (as 
estimated by IBirzan et~aH |2004| ). IBest et ail ((2006) esti- 
mated the time-averaged mechanical luminosity output as- 
sociated with radio source activity for each galaxy, as a func- 
tion of its mass. They found that, for elliptical galaxies of 
all masses, the time averaged radio heating almost exactly 
balanced the radiative cooling losses from the hot haloes of 
the ellipticals. They argued that the radio AGN feedback 
may therefore play a critical role in galaxy formation, pre- 
venting further cooling of gas onto elliptical galaxies and 
thus explaining why these galaxies are old and red (see also 
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Benson et al.l 200jj ; Scannapieco et al.l |2005l; ICroton et al.l 



2006; 



Bower et 



OOll; ._ ...... 

al.ll2006l ; ICattaneo et al.ll2006l ). 
iBest et aL l|2006T l also concluded that, unless brightest 
cluster galaxies showed a different mode of radio activity to 
that of ordinary elliptical galaxies, then they could not pro- 
vide sufficient heating to balance the cooling in clusters of 
galaxies. As discussed above, BCGs are found to be radio- 
luminous much more frequently than other cluster ellipti- 
cals, but the stron g mass depe n dence of the radio-loud AGN 
fraction found by Best et ail (I2005T ) means that this may 
solely be due to their very much higher masses. The goal of 
the current paper is therefore to investigate this issue, by de- 
termining both the radio-loud fraction, and the distribution 
of radio luminosities, for a large sample of brightest group 
and cluster galaxies with a wide range of masses and host 
group / cluster properties. This will then ascertain the im- 
portance of radio-AGN heating in clusters. A secondary goal 
of the paper is to investigate the role that non-BCG cluster 
galaxies may have in heating the intracluster medium. 

The layout of the paper is as follows. Section [2] describes 
the various data samples that are used for this analysis. An 
analysis of the radio properties of brightest group and clus- 
ter galaxies is carried out in Section^ while Section[3]inves- 
tigates the radio properties of the other group and cluster 
galaxies. These are compared against signatures of optical 
AGN activity in Section [S] The implications of these results 
for the origin and fuelling of low luminosity radio source ac- 
tivity are discussed in Section]!)] In Section[7]the radio-AGN 
heating rates are estimated in clusters and groups across 
a wide range of masses, and compared these with the ra- 
diative cooling rates. The implications of these results for 
AGN heating of cooling flows are discussed in Section [H] 
and conclusions are drawn in Section [9] Throughout the pa- 
per, the cosmological parameters are assumed to have values 
of n m = 0.3, Q.A = 0.7, and H = 70 kms _1 Mpc _1 . 



2 SAMPLE SELECTION 

2.1 Definition of the galaxy samples 

The Sloan D i gital Sky Survey <|York et all |2000| ; 
IStoughton et all l2002h is a five-band photometric and 
spectroscopic survey which will ultimately cover about 
a quarter of the extragalactic sky; the fourth data re- 
lease of this survey (DR4; lAdelman-McCarthv et all 
2006) includes spectroscopy of over half a million ob- 
jects. The parent sample for the current study is dr awn 
from the 'main galaxy sample' (jStrauss et al.l I2002T I of 
the SDSS DR4, comprising those galaxies with magni- 
tudes in the range 14.5 < r < 17.77. The spectra of 
these galaxies have been used to derive a large num- 
ber of physical properties, with catalogues of measured 
and derived parameters being publically available (see 
htt p:/ /www. mp a-garching.mpg.de/SDSS/ ). 

iMiller et all (120051 ) have used the SDSS data to identify 
clusters of galaxies as overdensities in a seven-dimensional 
space of position and colour. This 'C4' cluster catalogue is 
currently available for DR3, and identifies 1106 clusters with 
redshifts 0.02 < z < 0.16. Miller et al. also identify two can- 
didate brightest cluster galaxies for each of their clusters: 
the 'mean galaxy' is that closest to the peak of the galaxy 



density field, and the 'brightest galaxy' is the brightest spec- 
troscopically confirmed cluster member with a projected po- 
sition within 500/i _1 kpc (where h = H /1QQ kms _1 Mpc _1 ) 
of that density peak and a redshift within four times the ve- 
locity dispersion of the mean cluster redshift. In many cases, 
however, neither of these tu rns out to be the true BCG (see 
Ivon der Linden et al.l [20061 ). For example, in about 30% of 
dense clusters, the BCG is not included in the spectroscopic 
data, due to the problems of fibre collisions, and so is missed 
by the C4 catalogue. 

In a companion paper, Ivon der Linden et alj (2006) 
have re-analysed the C4 cluster catalogue in order to pro- 
vide a robust sample of clusters with well-defined BCGs. 
The reader is referred to that paper for a detailed discus- 
sion of how the cluster sample and BCGs were defined. In 
brief, the sample was first restricted to the 833 C4 clusters 
with redshifts z < 0.1. The brightest cluster galaxies were 
then identified using all available information (magnitude, 
colour, morphology, redshift if in the spectroscopic sample), 
together with visual inspection of colour images of the clus- 
ters (see von der Linden et al. for details). During this pro- 
cess, clusters which were identified to be substructures of 
larger clusters were removed from the sample, as were those 
for which an iterative algorithm to determine the cluster red- 
shift, velocity dispersion (ov) and virial radius (r2oo) either 
did not converge or retained three or fewer galaxies within 
3(J V and lr2oo of the cluster centre. This resulted in a final 
sample of 625 systems, for which the size, velocity disper- 
sion, BCG and cluster membership were all well determined. 
It turns out that some of these systems have velocity disper- 
sions which indicate that they are galaxy groups rather than 
clusters; the terms 'cluster' and 'brightest cluster galaxy' 
are used loosely in this paper, to include both clusters and 
groups. 

A further crit ical re-analysis carried out by 

Ivon der Linden et al. I l|2006h has been to correct the 
SDSS photometry for these BCGs. The SDSS photometry 
systematically underestimates the luminosities of nearby 
large galaxies, particularly in cluster environ ments, because 
it overestimates the leve l of sky background llBernardi et all 
l2006t lLauer et al.l 120061 ; Ivon der Linden et all 120061 ). This 
then leads to underestimates of the galaxy mass, which 
would impact upon the current studies, von der Linden 
et al. derived a method to correct the photometry for this 
sky over-subtraction; their method uses both the 'local' 
and 'global' sky estimates provided by the SDSS pipeline, 
combined according to a weighted ratio of the luminosity of 
the galaxy to that of its neighbouring galaxies. They show 
that this method is successful in reproducing magnitudes 
which match with independent measurements (see their 
paper for more details). 

In the current paper, the 'brightest cluster galaxy' sam- 
ple generally corresponds to the 625 BCGs defined from this 
sample of clusters. In 141 of these clusters, however, the 
BCG is not contained within the spectroscopic catalogue; 
therefore, where comparison is being made with properties 
derived from spectroscopic observations (e.g. emission line 
strengths), the sample is restricted to only those 484 BCGs 
with spectroscopic data. Since the BCGs with only photo- 
metric data have generally been excluded from the spectro- 
scopic sample due to random effects (especially fibre colli- 
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sions) it is not expected that this will produce any significant 
bias in the sample. 

The results for the BCGs are compared with those of 
'all galaxies'. The 'all galaxy' sample corresponds to those 
galaxies within the main galaxy sample of the SDSS DR4 
release, with redshifts 0.02 < z < 0.1; these redshift cuts 
match those applied to the cluster catalogue, and von der 
Linden et al. have also applied their photometric correction 
method to these galaxies in an equivalent way. The all galaxy 
sample is discussed by von der Linden et al.: it includes a 
mix of galaxy properties, in a wide range of environments, 
but in the luminosity range where the sample overlaps with 
that of the BCGs it is dominated by early-type galaxies. 

A 'non-BCG cluster galaxy' sample is also constructed. 
This consists of those galaxies within the 'all galaxy' sample 
which are within 3r2oo and 3<r v of one of the 625 clusters 
identified by von der Linden et al., but which are not iden- 
tified as a BCG. 

For all galaxies, the rest-frame luminosities were calcu- 
lated from the re -calibrated photomet ry using the KCOR- 
rect algorithm dBlanton et al. 2003), which determines 
the best composite fit to the observed galaxy fluxes of 
a large number of template stellar population spectra, 
including three alternative model s for dust extinction. 
Following lBlanton~fc Roweisl (|2007h . this same algorithm 
was used to derive stellar mass estimates: as shown by 
those authors, the stellar masses derived agree to typ- 
ically within 0.1 dex of those der i ved by other meth- 
ods (e.g. those of iKauffmann et alJ l2003ai , using the z- 
band magnitude and a mass-to-light ratio determined 
by the strengths of the 4000A break and the H6 ab- 
sorption in the galaxy spectrum). Other physical proper- 
ties of the galaxies, such as their emission line strengths, 
were taken from the MPA pipeline catalogues available at 
http:/ /www. mpa-garching.mpg.de/SDSS/ 

2.2 Identification of radio-loud AGN 



iBest et al. (2005) identified the radio-emitting galaxies 
within the main spectroscopic sample of the SDSS DR2, 
by cross-comparing these galaxies with the National Radio 
Astronomy Observa tory (NRAO) Very Large Array (VLA) 
Sky Survey (NVSS: ICondon et al.|[l99§ ) and the Faint Im- 
ages of_the_lljidto Skyat Twenty centimetres (FIRST) sur- 
vey l|Becker et al.|[l995T ). The use of a combination of these 
two radio surveys allowed a radio sample to be derived which 
was both reasonably complete (« 95%) and highly reliable 
(~ 99%). They then used the optical properties of the galax- 
ies to separate the radio-loud AGN from the radio-detected 
star-forming galaxies. This work has now been extended to 
include the DR4 data (Best et al. in prep), and these re- 
sults were used to determine which galaxies from the SDSS 
spectroscopic samples are radio-loud AGN. 

For the 141 brightest cluster galaxies with only photo- 
metric data, the same cross-comparison of the galaxy lo- 
cations with the NVSS and FIRST radio surveys was re- 
peated in order to determine which had associated radio 
emission. However, the lack of spectroscopic data for these 
galaxies meant that the previous method of distinguishing 
radio-loud AGN from star-forming galaxies could not be 
used. Instead, for the 37 photometric BCGs with radio de- 
tections, the classification was based solely upon the radio 
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Figure 1. The histogram shows the distribution of radio lu- 
minosities of the radio— detected BCGs with SDSS spectroscopic 
data, split by their classification as either radio-loud AGN (un- 
shaded) or radio-detected star— forming galaxies (shaded). The 
vertical lines at the top of the plot show the radio luminosities of 
the BCGs with only photometric data. 



luminosity. Figure [T] shows the distribution of radio lumi- 
nosities for the BCGs with spectroscopic data, split by their 
classification as either radio-loud AGN or radio-detected 
star-forming galaxies. Above I/i.4GHz = lO^WHz- 1 all 
but one radio detected BCGs are AGN, whereas below that 
luminosity nearly half are star forming galaxies. The vertical 
lines at the top of the plot indicate the radio luminosities of 
the BCGs without spectroscopic data. The 32 photometric 
BCGs with L 1.4GHz > lO^^WHz" 1 are almost certainly 
radio-loud AGN, and henceforth are classified as such. The 
other 5 all have radio luminosities below 10 22,6 WHz _1 and 
are likely to include a mixture; their classification is less 
critical because most analyses are limited to only the BCGs 
brighter than 10 23 WHz _1 , but conservatively it is assumed 
that these are all star forming galaxies. It should be empha- 
sised that the number of BCGs for which the origin of the 
radio emission is ambiguous is far too small for any misclas- 
sification to influence the results of the paper. 



3 THE RADIO ACTIVITY OF BRIGHTEST 
CLUSTER GALAXIES 

Using t he radio sources defined from the SDSS DR2 cat- 
alogue, IBest et all l|2005l ) showed that the probability of 
a galaxy to be a radio-loud AGN was dependent primar- 
ily upon its mass, determined either as the stellar mass 
(/radio-ioud oc M, 15 ), o r the black hole mass (/ ra dio-ioud oc 
Mbh)- As discussed bv lBest et alJ (|2005l ). the difference be- 
tween the slopes of these two relations arises mainly because 
of the increasing fraction of disk-dominated galaxies, with 
small black holes, at stellar masses M» < 1O U M . These 
host fewer radio-loud AGN, and thus decrease the radio- 
loud fraction at low stellar masses, leading to a steeper de- 
pendence on stellar mass. 

Figure[5]shows the percentage of brightest cluster galax- 
ies that are radio-loud AGN, as a function of stellar mass. 
For comparison, the equivalent relation is also shown for 
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the all galaxy sample. The brightest cluster galaxies are 
clearly offset from the 'all galaxy' relation, and their re- 
lation shows a different slop43 : the fraction of BCGs that 
is radio-loud increases roughly as Ml' and is an order of 
magnitude higher than that of all galaxies for stellar masses 
below 1O 11 M0, but comparable to that of all galaxies at stel- 
lar masses above 5 x 10 11 Mq . The increasing similarity of the 
two radio-loud fractions with increasing mass is actually not 
surprising, since at masses above 10 Mq around a third 
of the 'all galaxy' population are classified as BCGs, and 
many of the remainder may be BCGs missed by the current 
sample. Note that from these data it is not possible to tell 
whether the increase in the radio-loud fraction of BCGs, rel- 
ative to that of other galaxies, is due to radio-AGN activity 
being triggered more frequently in BCGs, or the BCG radio 
sources being longer lived (e.g. because the higher pressure 
environment of the dense cluster environment slows down 
the source expansion). 

Black hole mass may be considered to be a more 
fundamental parameter t han stellar mass when consider- 
ing nuclear activity, and iBest et al.l (j2005h estimated this 
for the SDSS DR2 galaxies using the velocity dispersio n 
versus black hole mass relation of iTremaine et all (2002): 
log(M BH /M ) = 8.13 + 4.021og(o- t /200kms _1 ). However, 
there is some doubt as to whether this 'normal galaxy' re- 
lation is also applicable for BCGs, or whether thes e follow 
a different relation (e.g. Ivon der Linden et aT]|2006t ). In ad- 
dition, velocity dispersion measurements are only available 
for those BCGs with spectroscopic observations. In order to 
avoid any potential biases therefore, and also to increase the 
sample size, in this paper the analysis is presented only as 
a function of stellar mass. As a check, relations were con- 
structed as a function of black hole mass, assuming the ve- 
locity dispersion to black hole mass conversion to be the 
same as for normal galaxies; these are in broad agreement 
with those presented here for stellar mass. 

Figure |3 compares the radio-loud fraction versus stellar 
mass relation for BCGs in rich versus poor clusters. The 
probability for a BCG of given stellar mass to host a radio- 
loud AGN is independent of the velocity dispersion of the 
surrounding cluster. Because BCGs in the richest clusters 
tend to be more massive than those in poorer clusters, the 
fraction of BCGs which host radio-loud AGN still turns out 
to be higher in more massive clusters (marginally at least; 
21 ± 2% in the ov > 500km s _1 clusters in this sample, 
compared to 17 ± 2% for those with a v < 500km s" 1 ), but 
it is the mass of the BCG, rather than the properties of 
its surroundings, which controls the likelihood of it being 
radi o-loud. 

IBest et all (|2005h also found that for 'all galaxies' the 



It should be stressed that the same conclusions are reached if 
the plots are made as a function of galaxy luminosity, i nstea d 
of stellar mass, as demonstrated in Ivon der Linden et al. ] tod) . 
Similarly, if the 'all galaxy' sample is restricted to just those 
classified as early-type galaxies by using their SDSS morphologi- 
cal parameters concentr ation index and surface mass density (cf. 
iKauffmann et al]|2003bl , and discussion therein) then the results 
are also largely unchanged at masses above 5 X 10 10 Mq, where 
it overlaps with the BCG sample. This is because the galaxies at 
high stellar masses are dominated by old ellipticals, whose mass- 
to-light ratios do not vary strongly. 
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Figure 2. The percentage of galaxies that are radio-loud AGN, 
as a function of stellar mass, for 'all galaxies' and for 'brightest 
cluster galaxies'. Brightest cluster galaxies are more likely to be 
radio-loud than other galaxies of the same stellar mass, and the 
two relations have different slopes. 
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Figure 3. The percentage of brightest cluster galaxies that are 
radio-loud AGN, as a function of stellar mass, split by the veloc- 
ity dispersion of the clusters. At given stellar mass, the radio— loud 
fraction of BCGs is indistinguishable between rich and poor clus- 
ters. 



distribution of radio luminosities of the radio-loud AGN was 
independent of the mass of the host galaxy. This result is 
shown in the top panel of Figure [3] for three different mass 
ranges. Also shown on the same plot are the distributions of 
radio luminosities for brightest cluster galaxies in the same 
mass ranges. The shape of the radio luminosity function 
does not change significantly either as a function of mass for 
brightest cluster galaxies, or between the brightest cluster 
galaxy and all galaxy sample^]: only the normalisation of 
the radio luminosity function changes. The bottom panel of 



2 The only variation in shape is a flattening at low radio luminosi- 
ties, of both BCGs an d all galaxies, whe n the radio-loud fraction 
exceeds 20 - 30%; as IBest et al.l J2005l> argued, such flattening 
must occur at some point, since the radio-loud fraction cannot 
exceed 100%. 
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Figure 4. Top: The percentage of galaxies which are radio— loud 
AGN brighter than a given radio luminosity, for 'all galaxies' 
(black lines) and the subsample of 'brightest cluster galaxies' (red 
lines) in three different ranges of stellar mass. There is no evidence 
for any significant difference in the shape of the luminosity func- 
tions, either as a function of stellar mass, or between brightest 
cluster galaxies and all galaxies: only the normalisation of the 
relation changes. Bottom: The distribution of radio luminosities 
of BCGs split by the velocity dispersion of their host clusters, 
for BCGs with masses in the range 10 11 < M*/M® < 5 x 10 11 , 
compared to that of normal galaxies. 



Figure [3] compares the distribution of radio luminosities of 
BCGs in rich and poor clusters: again, these are broadly in 
agreement, although there is tentative evidence that there 
may be more high luminosity radio sources in richer clusters. 
A larger sample will be required to confirm this result. 



4 RADIO-LOUD AGN ACTIVITY IN OTHER 
CLUSTER GALAXIES 

Many authors have examined the environments of radio- 
loud AGN (e g.lPrestage fc Peacockllf98sl ; Iffill fc Lillv||l99ll : 
iMiller et ail |2002| ; iBestJ 12004! ) and have found that these 
sources appear to favour galaxy group and weak cluster 
environments. Nevertheless, Ledlow fc Owenl (| 1996f ) demon- 
strated that the bivariate radio-optical luminosity function 
of cluster galaxies is very similar to that of field galaxies. 
This indicates that the preference for radio sources to be 



Figure 5. The percentage of galaxies that are radio-loud AGN, 
as a function of stellar mass, for 'non-BCG cluster galaxies' (red 
line) as compared to the 'all galaxy' (dotted black line) and 
'brightest cluster galaxy' (solid black line) samples. The other 
coloured lines indicate the radio-loud AGN fractions for subsets 
of the non-BCG cluster galaxies within different ranges of pro- 
jected cluster-centric radii: the green, orange and cyan lines show 
the relations for galaxies with redshifts within 2<r v of the mean 
cluster rcdshift and with projected radii of 0— 0.2r200i 0.2-0.5r2oo 
and 0.5-1.0r2oo respectively, while the dark blue line shows the 
equivalent relation for cluster galaxies outwith 1.0r2oo- 



located in dense environments may simply be due to the 
concentration of massive galaxies within clusters, coupled 
with the strong mass dependence of the radio-loud fraction. 
The SDSS data and the C4 cluster sample allow this issue 
to be investigated in more detail. 

Figure[5]shows the fraction of non-BCG cluster galaxies 
(red line) that are radio-loud AGN, as a function of stellar 
mass, in comparison to all galaxies and to brightest cluster 
galaxies. The radio-loud AGN fraction of cluster galaxies 
matches that of all galaxies at all stellar masses (except pos- 
sibly the lowest mass bin, where it might be slightly lower). 
Figure [6] similarly shows that the distribution of radio lumi- 
nosities of the cluster galaxies is the same as that of non- 
cluster galaxies. Overall, then, for galaxies of given stellar 
mass the radio luminosity functions are the same ins i de an d 
outside of clusters, as was found bv lLedlow fc Owenl |l996). 
They are also independent of the velocity dispersion of the 
host cluster, as shown in Figure [7] 

Figures [5] and [6] also show the radio-loud AGN frac- 
tion and radio luminosity distribution for various subsets of 
cluster galaxies in different ranges of projected radius from 
the cluster centre. It is striking that galaxies with redshifts 
within 2<t v of the mean cluster redshift and 0.2r2oo of the 
cluster centre have an enhanced probability of being radio- 
loud AGN, approaching that of the BCGs, although once 
again with the same distribution of radio luminosities. Out- 
side of 0.2r2oo there is no substantial change in the radio- 
loud fraction with radius (there may be a marginal fall with 
increasing radius, but this is within the errors, and insignif- 
icant compared to the increase within 0.2r2oo). 

One concern is that the increase in the radio-loud frac- 
tion of galaxies near the centre of the cluster may arise as 
a result of mis-identification of some BCGs. The evidence 
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Figure 6. The percentage of galaxies which are radio— loud 
AGN brighter than a given radio luminosity for 'non— BCG clus- 
ter galaxies' (red line) as compared to the 'all galaxy' (dotted 
black line) and 'brightest cluster galaxy' (solid black line) sam- 
ples. These plots are constructed considering all galaxies in the 
10 11 < M»/Af < 10 12 mass range, but then scaling the BCG 
line uniformly down by a factor of 1.7 to account for the higher 
median mass of BCGs within this mass range (the median masses 
of galaxies in all other subsamples are statistically indistinguish- 
able). The other coloured lines indicate the equivalent relations 
for different subsets of the non-BCG cluster galaxies with differ- 
ent projected radii, as defined in Figure [5] The radio-loud AGN 
fraction is boosted for galaxies within 0.2r2oo of the cluster centre, 
but retains the same overall shape. 
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Figure 7. The distribution of radio luminosities of non-BCG 
cluster galaxies, split by the velocity dispersion of their host clus- 
ters, as compared to all galaxies and to brightest cluster galaxies. 
Only galaxies with masses in the range 10 11 < M„/Mq < 5 x 10 11 
are included on the plot. 



suggests, however, that this is not the case. There are 13 
clusters where a radio source is associated with a non-BCG 
cluster galaxy within 0.2r2oo- Of these, the identified BCG 
is also radio-loud in 6 (= 46±19%) of the clusters: this frac- 
tion is comparable to (or even above) that of other BCGs 
of the same mass (~ lO n ' 3 M0), and far above that of non- 
BCG cluster galaxies. Further, of the 7 clusters with a radio- 
loud galaxy within 0.2r2oo but a radio-quiet BCG, visual 



analysis indicates that the BCG identification is absolutely 
unambiguous in 5 of the cases, and should be reliable in 
the other two. Finally, the distribution of ratios of radio 
galaxy mass to BCG mass is indistinguishable between the 
13 r < 0.2r2oo radio galaxies and the radio galaxies at larger 
radii, in contrast to what would be expected if these were 
misidentified BCGs. In conclusion, therefore, the enhanced 
radio-loud fraction for galaxies within 0.2r2oo of the cluster 
centre appears to be a genuine physical effect. 



5 EMISSION-LINE AGN ACTIVITY IN 
CLUSTER GALAXIES 

It is interesting to compare the enhanced radio-loud AGN 
activity of brightest cluster galaxies and other galaxies in 
the centra l 0.2r2oo of the cluster, with their emission line 
properties. iKauffmann et al. (2003c) identified emission-line 
galaxies in the SDSS spectroscopic sample, and used the 
[OIII] 5007 / H/3 vers us [Nil] 6583 / Ha e mission line ra- 
tio diagnostic diagram l|Baldwin et al.lll98ll) to separate out 
galaxies with AGN activity from those galaxies where the 
emission lines are associated with star formatiorJf|. The frac- 
tion of galaxies hosting emission line AGN rises shallowly up 
to 1O 1O ' 5 M0, and is then flat at higher masses, which is sub- 
stantially different from the v ery steep mass d ependence of 
the radio-loud AGN fraction |Best et al.ll2005l ). This result 
is shown in Figure |S]for t he 'all galaxy' sample (fo r AGN de- 
fined in the same way as IKauffmann et al.ll2003d . and with 
£[oiii]5007 ^ W 6 ' 5 Lq) and is compared to the fraction of 
galaxies hosting emission-line AGN for the brightest cluster 
galaxy and non-BCG cluster galaxy samples. Also shown are 
the relations for subsamples of cluster galaxies at different 
cluster-centric radii. Note that once again, making this plot 
as a function of galaxy luminosity instead of mass leads to 
essentially the same results. 

Brightest cluster galaxies are found to be less likely to 
possess emission-line AGN activity than 'all galaxies' of the 
same stellar mass, by a factor of 2-3 at all masses. This 
suppression of emission-line AGN activity is also true of 
other cluster galaxies: cluster galaxies as a whole are found 
to be about 20% less likely to show emission-line AGN ac- 
tivity than all g alaxies, consistent with the results of prev i- 
ous studies (e.g. iMiller et al.ll2003l : IKauffmann et alj|2004l ). 
Splitting the cluster galaxies into subsamples at different 
radii demonstrates that the suppression gets progressively 
stronger as the cluster centre is approached: outside T2oo 
the fraction of cluster galaxies displaying emission-line AGN 
activity matches that of 'all galaxies', but this drops by a 
factor of at least 2 for those galaxies within 0.3r2oo- Inter- 
estingly the BCGs show the same likelihood of emission line 
AGN activity as other galaxies of the same mass towards the 
centre of the cluster: unlike for the radio activity, the spe- 
cial location of BCGs does not appear to lead to different 
emission-line AGN activity. 

This result may seem to be at variance with the lu- 
minous extended emission line structures known to exist 



3 It should be stressed that since the parent sample for that anal- 
ysis, like the current one, was the SDSS 'main galaxy sample', the 
AGN selected are mostly type-II AGN. Powerful type-I quasars 
will be missed from the sample. 
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around some BCGs (e.g. iHeckman et al.lll989l ; lHatch et alj 
120061 . and references therein). Part of this difference may be 
because the SDSS fibre apertures are small (~ 5kpc diam- 
eter at 2 ~ 0.1), whilst the line emission seen in the above 
studies can extend for several tens of kiloparsec. More im- 
portant, however, is a difference in the type of BCG studied. 
lEdwards et~ai] (|2007T l recently studied the line emission from 
BCGs and found that the fraction of the BCG population 
as a whole that display emission lines is about 15%, and is 
comparable to that of other massive galaxies near the clus- 
ter centres. These results are in full agreement with those 
found here, shown in Figure [8] However, these authors also 
find that if they restrict analysis to just those BCGs which 
are within 50 kpc of the X-ray centre of a cooling flow clus- 
ter, then the fraction which display emission lines rises to 
~ 75%. ft is such massive, strong cooling flow clusters that 
have been the main focus of the previous emission line stud- 
ies, and the extended emission line activity detected is asso- 
ciated with the cooling flow. The current results and those 
of Edwards et al. demonstrate that more typical BCGs do 
not have enhanced emission line activity. 

A further interesting test concerns the relat ionship be- 
tween radio and optical emission line activity. iBest et all 
(2005) found that for 'all galaxies' the probability of a galaxy 
hosting a low luminosity radio source was independent of 
whether or not it was an emission-line AGN. Figure [9] in- 
vestigates this for the brightest cluster galaxies and the non- 
BCG cluster galaxies. Unlike the situation for all galaxies, 
the radio-loud fraction for brightest cluster galaxies does de- 
pend upon the emission line properties: BCGs which show 
emission-line AGN activity are 2-3 times more likely at all 
masses to host radio-loud AGN than those which are op- 
tically inactive — although the latter still show a higher 
radio-loud fraction than non-BCG galaxies. Other cluster 
galaxies behave in the same way as 'all galaxies', with the 
radio-loud fraction being the same regardless of the optical 
properties. 



6 IMPLICATIONS FOR THE ORIGIN OF LOW 
LUMINOSITY RADIO-LOUD AGN 

The apparent independence of emission line and radio- 
loud AGN activity, coupled with the fundamentally differ- 
ent dependencies of their prevalences on galaxy mass, led 
IBest et al.l l|2005h to conclude that low luminosity radio- 
loud AGN activitjQ and emission line AGN activity are dis- 
tinct physical phenomenon. They argued that emission-line 
AGN activity tracks the accretion of cold material onto black 
holes through standard accretion disks, and thus reflects the 
major mode of growt h of the black holes (see discussion in 
IHeckman et al.ll2004 ). whilst low luminosity radio sources 
represent re-triggering of pre-formed massive black holes 
in old elliptical galaxies. The accretion rates in these ra- 
dio sources are comparably low, and the activity is reflected 



4 The volume probed by the SDSS data is insufficient to investi - 
gate the rare, powerful, classical double jFanaroff & Riley dl974h 
class 2, or FRII) radio sources, which do have associated optical 
AGN activity, but instead only investigate the distinct population 
of lower luminosity (mostly FRI) radio sources. 
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The percentage 

'[OIII]5007 



of galaxies that show emission-line 
AGN activity (with £[oin]5007 ^ 10 6 ' 5 ij©), as a function of stel- 
lar mass, for 'all galaxies' (dotted black line), 'brightest cluster 
galaxies' (solid black line) and 'non-BCG cluster galaxies' (red 
line). Also shown are equivalent relations for subsamples of clus- 
ter galaxies in different cluster-centric radial ranges. Emission line 
AGN activity is suppressed within r200 of the cluster centre, with 
the strength of the suppression increasing with decreasing radius. 
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Figure 9. The percentage of brightest cluster galaxies which 
are also emission— line selected AGN (solid red line) and those 
which are optically— inactive (dashed red line), that are classified 
as radio-loud AGN, as a function of stellar mass. This is compared 
to the fraction of normal galaxies that are radio— loud AGN (black 
line; this is independent of whether or not they optical AGN), and 
with the equivalent relations for non-BCG cluster galaxies (blue 
lines). 



in the production of low power radio jets but relatively lit- 
tle emission at op tical, ultraviolet and infrared wavelengths. 
IBest et alJ <|2005f ) argued that the accreting material in the 
radio mode was hot gas, cooling out of the envelopes of el- 
liptical galaxies. The results for the brightest cluster galax- 
ies lend further credence to this scenario: since BCGs are 
located at the bottom of the cluster potential well, they 
are also the repository for cooling intracluster gas. This in- 
creased supply of hot gas naturally leads to the observed 
increased radio-loud AGN fraction. 
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iBest et all (|2006l ) suggested that one possibility for the 
accretion mechanism is Bondi-Hoyle accr etion from a stron g 
cooling flow (see also the recent work of I Allen et al,l [2006): 
the accretion rate for this is given by ICroton et al.l (2006) 
as mBondi ~ Gfj,m p kTMBn/A, where fxm p is the mean mass 
of particles in the gas, and A is the cooling function. For 
isothermal elliptical galaxy halos, T <x a 2 ; the black hole 
mass versus velocity dispersion relation fe.g. lTremaine et al.l 
2002) then implies T oc . A is relatively independent of 
temperature (and hence black hole mass) at the temperature 
of elliptical galaxy haloes, and therefore the Bondi accre- 
tion rate scales roughly as M^, which is approximately the 
same exponent as the radio-loud AGN fraction for all galax- 
ies. Interestingly, for BCGs the approximation T oc Af£^ is 
not valid, because the temperature of the gas accreting on 
to the BCG is not directly related to the properties of the 
host galaxy, but rather is controlled by the larger-scale en- 
vironment. In a cluster of given temperature, therefore, the 
Bondi-Hoyle accretion rate might be expected to scale lin- 
early with the black hole mass (and hence stellar mass) of 
the BCG, which is exactly what is observed in Figure [5] 

One surprising result is that there is no increase in 
the radio-loud AGN fraction in higher velocity dispersion 
(mass) clusters: these would be expected to have higher cool- 
ing rates, raising the radio-loud fraction. For the Bondi- 
Hoyle mechanism, at fixed black hole mass the accretion 
rate scales as T/A, and at the temperature of clusters, A 
increases approximately as T 0,5 . The accretion rate there- 
fore scales as T ' 5 , and thus roughly linearly with the cluster 
velocity dispersion, a v . This would predict nearly a factor 
of two difference between the accretion rates, and hence the 
radio-loud AGN fractions, of the clusters in the low and 
high velocity dispersion bins, which ought to be (just) de- 
tectable from the current data. The failure to observe this 
effect may be because the relevant temperature is that of 
the gas in the central regions of the cluster, and most mas- 
sive clusters have cooler cores associated with the cooling 
flow; this reduces the expected difference in accretion rate 
between high and low mass clusters, possibly to within the 
errors of the current measurements. Investigating this with 
a much larger sample of clusters, or using clusters for which 
the core temperatures of the X-ray gas have been measured, 
would provide a critical test of this model. 

The increased likelihood of radio-loud AGN activity 
for non-BCG galaxies within 0.2r2oo, particularly the more 
massive ones, suggests that the accretion of gas onto these 
galaxies is also controlled to some extent by the intracluster 
medium. It is intriguing that the cooling radius of clusters 
(that is, the radius within which the cooling time is less than 
the Hubble time) is typically about 10% of r-2oo- thus, the 
intracluster medium affects the radio-loud fraction of only 
those galaxies within, or close to, the cooling radius. 

The emission-line properties of the BCGs provide fur- 
ther support for the interpretation that different accreting 
material is responsible for low-luminosity radio-loud AGN 
activity than that for emission line AGN activity. In contrast 
to the radio-loud AGN activity, the emission-line properties 
of BCGs are in general suppressed relative to galaxies out- 
side clusters, and are comparable to those of other galaxies 
of the same mass in the inner regions of the clusters. This 
implies that there is a lack of cold gas to fuel the black 
holes in these galaxies. This result is undoubtedly related 
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Figure 10. An emission line diagnostic plot for the brightest 
cluster galaxies with emission lines, divided into different classifi- 
cations for the source of ionisation of the emission lines following 
iKauffmann et alj<2003ch . The majority of these BCGs, especially 
the radio-loud ones, have LINER-like emission line spectra. 

to the decreased leve ls of star formation activity in clus- 
ter environments (e.g.lLewis et al ] |2002l : Gome z et al 1 l2003l : 
iMateus fc Sodre]|2004l . and references therein). 

Interestingly, those BCGs which do show emission-line 
activity also show an enhanced level of radio-loud AGN ac- 
tivity, suggesting that in BCGs (in contrast to other galax- 
ies) there may be some connection between the processes 
which drive the radio and emission-line activity. This is 
consistent with the discussion in Section [5] that the BCGs 
of cooling flow clusters frequently have associated extended 
emission-line nebulae. The properties of these are simi- 
lar to those of low-ionisation nuclear emission-line regions 
(LINERs; e.g. IVoit fc Donahulll997r ). and so may be mis- 
taken for AGN activity. Thus, the emission lines detected 
in some BCGs may not actually arise from AGN activity, 
but rather from the cluster-scale cooling flow which fuels 
the radio-loud AGN activity. Figure [TO] shows the location 
of the SPSS emissio n-line BCGs on the traditional 'BPT' 
l|Baldwin et al.lll98ll ) emission line diagnostic diagram. As 
a guide, this diagram has been di vided into regions (as de- 
fined by IKauffmann et al.l |2003cJ) indicating the expected 
locations of star forming galaxies, Seyferts, LINERs, and 
galaxies whose spectrum is a star formation - AGN com- 
posite. The majority of the emission-line BCGs, especially 
those which are radio-loud AGN, are found in the LINER 
region of the diagram, and so some of this line emission may 
indeed be directly associated with the cooling flows. 



7 RADIO-LOUD AGN HEATING OF 
COOLING FLOWS 

7.1 The mass— dependent heating rate of BCGs 

Radio-loud AGN have relatively short lifetimes (10 7 — 10 8 
years), and so the very high fraction of galaxies which host 
radio-loud AGN (over 30% at the highest masses) implies 
that this activity must be constantly re-triggered. The prob- 
ability of a galaxy of a given mass hosting a radio source of 
given radio luminosity (cf. Figure [4} can therefore be inter- 
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preted, probabilistically, as the fraction of its time that a 
galaxy of a given mass spends as a radio source of given lu- 
minosity. Monochromatic radio luminosity represents only a 
tiny fraction of the energetic output of a radio source, how- 
ever, with the mechanical energy of the radio jets being 2-3 
orders of magnitude larger. In order to determine the time- 
averaged heating output of radio-loud AGN, therefore, it is 
necessary to derive a conversion between the radio luminos- 
ity and mechanical luminosity of the radio sources. 

Any such conversion is necessarily approximate because 
the precise ratio between mechanical and radio luminosities 
is both unknown, and varies throughout the lifetime of a 
radio source. Nevertheless, it is reasonable to assume that 
there may be some statistical average relationship (e.g. aver- 
aging over all ages of a radio outburst), such that a conver- 
sion between the radio and mechanical luminosity functions 
can still be carried out at a population level. Such a conver- 
sion should be largely independent of the host galaxy mass 
(indeed, if it were not, then the similarity of the radio lu- 
minosity shapes for galaxies of different masses would be a 
remarkable co-incidence). The mean conversion relation can 
be determined empirically using data from the observed bub- 
bles and cavities produced by radio sources in the intraclus- 
ter me dium of groups and clusters of galaxies. iBirzan et al.l 
l|2004 ) compiled data for the energies, ag e s and radio lumi- 
nosities of these cavities, and iBest et all (|2006T > used these 
data to determine that: 

r / T \ 0.40±0.13 

Given the result of Section [3] that the shape of the ra- 
dio luminosity function is the same between galaxies of all 
masses, then converting the radio luminosity function to a 
mechanical luminosity function using Equation[TJmeans that 
the shape of the mechanical luminosity function is also sim- 
ilar for galaxies of different masses (and for BCGs). Inte- 
grating across this, Best et al. estimated the time-averaged 
AGN heating output of all galaxies, as a function of black 
hole masijf]: 

B M = 1.6 x 10 34 (Mbh/10 8 M s ) 1 - 6 W. (2) 

In this equation, the black hole mass dependence arises 
from directly from the scaling of the radio-loud fraction 
with black hole mass (ie. that more massive black holes are 
switched on for a larger fraction of the time). Errors in the 
conversion of radio to mechanical luminosity, and other such 
uncertainties, only affect the normalisation factor. This can 
be accounted for by introducing a factor / into Equation [5] 
giving: 

B M = 1.6 x 10 M /(M B H/10 8 Af©) 1 - 8 W. (3) 

The factor / accounts for a range of possible uncertain- 
ties, and these are discussed in detail in Section 17.21 It is 
likely, however, that i t has a valu e reaso nably close to unity. 
It is noteworthy that IBest et all (2006) found that a value 



/ ~ 1 was appropriate for heating to balance radiative cool- 
ing losses in elliptical galaxies. 

The method described above can also be applied to the 
brightest cluster galaxies. Taking the result from Section [3] 
that the radio luminosity function of BCGs has the same 
shape as that of 'all galaxies[f] and that only the normalisa- 
tion differs, then integrating the distributions in Figure |4l 
and incorporating the same factor / gives: 

JTbcg = 2.3 x lO 35 /(M,/lO n M )W. (4) 

As shown in Section [4] the time averaged heating rate 
of the non-BCG cluster galaxies outwith 0.2r2oo is the same 
as that for 'all galaxies', given in Equation^ For the cluster 
galaxies within 0.2r2oo, the boosted radio-loud AGN frac- 
tion will produce a heating rate between that of the BCGs 
and that of all galaxies, but the sample is too small to ac- 
curately determine this. For the analysis in the following 
sections, the limiting cases of H a \\ and Hbcg were both con- 
sidered for these galaxies, and it was found that there was 
negligible difference between the results produced. 

7.2 The uncertainty factor / 

The analysis method adopted here assumes that a good 
mean conversion can be made between radio and me- 
chanical luminosity. There are, however, examples of 
cluster radio sources ( e.g., MS0735.6+7421 Abell 1835; 
iMcNamara et all 120051 : IMcNamara et all 120061 ) where the 
mechanical power estimate is 2-3 or ders of magni t ude h igher 
than that estimated by Equation l4l iBirzan et al.l (|2004) and 
iRaffertv et ail (|2006h showed that the mechanical luminosity 
also relates strongly to the X-ray luminosity of the system, 
suggesting that a full understanding of these systems will 
require consideration of more parameters than simply the 
radio luminosity. Nevertheless, the uncertainty factor / in 
Equations [3] and [4] can be used to broadly encompass these 
additional effects. It is therefore worthwhile considering in 
detail what uncertainties are contained within the factor /, 
as this has important consequences for the interpretation 
of the results. The uncertainty can be categorised in three 
parts. 

(i) Uncertainty in the estimates of cavity mechanical lu- 
minosities. The calculations of Birzan et al. have significant 
uncertainties in the estimates of both the energies and the 
ages of the radio source cavities. In addition, Birzan et al. 
calculate the pV energy of the cavity, but the enthalpy of 
the cavity is given by ;^rj-pV, which is 4pV for the relativis- 
tic plasma of the radio lobes. There may also be additi onal 
heating directly from the radio jets: iNusser et all (|2006l ) ar- 
gue that the mechanical energy may exceed pV by a factor 
of 10. 

(ii) Uncertainty in the radio to mechanical luminosity 
conversion. As discussed above, the varying radio lumi- 
nosity during a source lifetime means that this conversion 



5 It is worth noting that iNipoti fc Binnevl ||2005| 1 adopt a simi- 
lar method but with the opposite approach: beginning from the 
assump tion that AGN mech anical heating balances cooling, they 
use the IBirzan et al. results to derive a radio luminosity 

function for radio-loud AGN, which they show to be in good 
agreement with that observed. 



6 Note that if the radio luminosity function for BCGs shows a 
weaker break (which would be consistent with the observations 
for clusters with tr v > 500km s —1 ) then the derived heating rate 
would be higher — but even if there was no break at all out to an 
extreme luminosity of 10 30 WHz _1 the increase would only be a 
factor of 1.7. 
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is only accurate in an average sense over the whole radio 
source population; however, to first order it ought to be 
both reasonable, and largely independent of host galaxy 
mass. There have been arguments that this conversion might 
be dependent upon environment, however: for fixed jet ki- 
netic power, radio luminosities could be higher in clus- 
ters, due to the confining effect of the dense intracluster 
medium reducing adiabatic expansion losses in the radio 
lobes and therefore increasing the radio synchrotron emis- 
sion (e.g. Bart hel fc Arnaudl 19961 . and references therein). 
If correct, this would lead to slightly lower values of / in 
denser systems. Alternatively, as suggested by systems such 
as MS0735. 6+7421 and Abell 1835 discussed above, individ- 
ual outbursts with extreme mechanical to radio luminosity 
ratios may become more important in clusters, effectively 
raising the value of /. 

(iii) Uncertainty in how much of the cavity mechanical 
energy gets converted to heat within the cooling radius. The 
heating rate calculations above assume that all of the me- 
chanical energy of the cavities is used to heat the intracluster 
medium, and that such heating occurs within the cooling 
radius of the cluster. If either the radio source energy is 
not efficiently converted to heat (e.g. because of inefficient 
production of sound/shock waves), or much of that heating 
occurs beyond the cooling radius of the cluster, then this 
could lead to values of / <C 1. 



7.3 The importance of the BCG in radio— mode 
heating 

The top panel of Figure [TT] shows the ratio of the time- 
averaged radio-mode heating rate, summed across all non- 
BCG cluster galaxiefl, to that of the BCG, for all of the 
625 SDSS clusters. These are calculated using Equations 
and 21 assuming that the same value of / is appropriate for 
both. In most groups and low mass clusters the heating is 
dominated by the BCG alone, but for the richer clusters, the 
sum total of other galaxies can provide significantly more 
heating than the BCG alone. This heating from the non- 
BCG galaxies is, however, spread over a volume very much 
larg er than the c ooling radius of the cluster. 

iPeres et alj (| 19981 ) have determined the cooling radii of 
an X-ray flux limited sample of 55 nearby clusters, and these 
cooling radii (converted to the cosmology adopted here) are 
plotted as a function of cluster velocity dispersion in the 
upper panel of Figure 1121 As most of these clusters are of 
high velocity dispersion, these data have been supplemented 
by three lower velocity dispersion groups for which coolin g 
radii ha ve also been mea sured: NGC155 dSun et al.ll2003h. 
RGH 80 (|Xue et ahlboolh and NGC6482 (jKhosroshahi et al l 
120041 '). The solid line shows a linear fit (in log space) to the 
data, calculated using the EM and Buckley-James linear re- 
gression techniques for censored data, within the ASURV 



7 Although the SDSS data only include the more luminous clus- 
ter galaxies, the very strong mass dependence of the radio— loud 
fraction means that the missing fainter galaxies should have a 
negligible contribution to the radio AGN heating rate. The SDSS 
data will also miss some of the brighter non-BCG cluster galax- 
ies, for example due to fibre collisions, but the number of these 
will be small. 
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Figure 11. Top: the ratio of the time-averaged radio-mode heat- 
ing rate of all non-BCG cluster galaxies, to that of the BCG, as 
a function of velocity dispersion, for the 625 SDSS clusters. Bot- 
tom: the same plot, but restricted to only those non-BCG cluster 
galaxies which are projected to lie within the cooling radius of 
the cluster. 



survival analysis package (jLaVallev et alj Il992l i: r CO oi = 
67(o"„/500kms -1 )°' 45 kpc. Using this fit to estimate the cool- 
ing radius of each SDSS cluster, the lower panel of Figure [TT1 
shows the ratio of radio-mode heating from non-BCG clus- 
ter galaxies projected within the cooling radius to that of 
the BCG. This analysis demonstrates that only the BCG is 
really relevant when calculating radio-mode heating within 
the cooling radius: unless there are other effects not consid- 
ered here (e.g. the efficiency factor / varies between BCGs 
and non-BCGs for some unknown reason) then the model of 
iNusser et al. I (|2006l ) whereby cooling flows are suppressed by 
the combined activity of all cluster galaxies does not seem to 
work. This is a result of the very strong mass dependence of 
the radio-loud AGN fraction (with BCGs being the highest 
mass galaxies) coupled with the additional boosting of the 
radio-loud fraction amongst the BCG population. 



7.4 AGN heating versus cooling in clusters 

Figure [TS] shows the time-averaged mechanical luminosity 
associated with radio-loud AGN, assuming f = 1, summed 
across all cluster galaxies, as a function of cluster veloc- 
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Figure 12. Top: the obse rved cooling radii o f clusters and group s 
of galaxie s , tak en fr om iPeres et al 1 Jl99ST). ISun et al.l hoo^t) , 
IXue et alj (12004!) and iKhosroshahi et al. I (l2004h~ as a function of 
cluster velocity dispersion. The solid line represents a fit to the 
data using survival analysis techniques to properly account for 
the censored data. Bottom: the fraction of the X— ray luminosity 
of those clus ters that arises fro m within the cooling radius, as 
calculated bv lPeres et al. (1998). 



ity dispersion. For comparison, four recently derived rela- 
tions between bolometr ic X-ray luminosity and velocity dis- 
persion are also shown (IXue fc Wull2000l; IMahdavi fc Gellerl 



l200ll ; lOrtiz-Gil et~aT]|2004l ; IPopesso et alj|2005h . The'com- 
bined heating rate of all radio-loud AGN increases much 
more slowly with cluster velocity dispersion than the a 4 
dependence of the radiative cooling rate. It is clear that 
(for / = 1) radio-loud AGN heating falls short of balancing 
the radiated X-ray luminosity for essentially all clusters and 
groups with a v > 300km s _1 , and by an order of magnitude 
at the high mass end. However, much of the X-ray luminos- 
ity of the clusters arises from radii larger than the cooling 
radius, and does not need to be balanced by heating in or- 
der to avoid catastrophic cooling collapse. It is instructive, 
therefore, to calculate what fraction of the X-ray luminosity 
within the cooling radius the BCG radio-mode heating is 
able to balance. 

IPeres et al.l (1 19981 ) determined the fraction of the to- 
tal X-ray luminosity that is emitted from within the cool- 
ing radius for their sample of 55 clusters. These fractions 



Figure 13. The data points show the time-averaged heating 
rate associated with radio-loud AGN activity, assuming / = 1 
in Equations \3\ and [4] summed across all galaxies in the clusters, 
as a function of cluster velocity dispersion. The solid line through 
these points, labelled '/ = 1', is a straight-line fit to these points. 
The two dotted lines parallel to this represent scaled versions of 
this fit, showing where the data points would be located if / = 0.1 
or / = 10. The remaining 4 lines indicate observed relations be- 
tween bolometric X-ray lumi nosity and velocity di spersion for 
clu ster and groups: sol i d line — lOrtiz-Gil et al.l j2004h; dotted line 
— IMahdavi fe Gellerl feoplll; dashed line — IXue k, Wul fcOOOT) ; 
dash-dot line — IPopesso et al ] <j2005l l. 



are plotted as a function of velocity dispersion in the lower 
panel of Figure 1121 The mean fraction of X-ray luminos- 
ity that arises from within the cooling radius is ~25%, 
and there is no significant trend with velocity dispersion. 
There is clearly considerable scatter around this mean value, 
but this is comparable to the scatter in the calculated 
AGN heating rates shown in Figure 1131 and so adoption 
of the mean value is reasonable. Combining this 25% frac- 
tion with t he X-ray luminosity t o velocity dispersion re- 
lation from lOrtiz-Gil et all l|2004| j g l gives L x (r < r coo i) w 
1.3 x lO^av/lOOOkms-y^W. 

Figure [T3] shows the fraction of the radiative cooling 
losses within the cooling radius that can be balanced by 
AGN heating from galaxies within that volume, assuming 
f — 1. Also shown are the location that the data points 
would occupy for alternative values of /. One conclusion is 
immediately apparent from these results: for a single value 
of /, radio-mode heating cannot balance radiative cooling 
losses for all systems from groups to clusters. 



8 This is the middle-most of the four X-ray luminosity to velocity 
dispersion relations shown on Figure [13] these relations are fitted 
to a combination of (mostly) clusters and (a few) groups. There 
is no consensus in the literature as to how the slope of this rela- 
tion changes for groups, and so here this mean relation is simply 
extrapolated down to the group regime; the conclusions of the 
paper would be unaffected even if the slope changed significantly 
below 300km s velocity dispersion. 
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Figure 14. The data points show the ratio of the time— averaged 
heating rate associated with radio-loud AGN activity of all galax- 
ies within the cooling radius, assuming / = 1 in Equations [3] 
and[4] to the radiative cooling losses within the cooling radius, as 
a function of cluster velocity dispersion. The solid line through 
these points, labelled '/ = 1', represents the fit to these points. 
The dotted lines are scaled versions of this fit, showing where the 
data points would be located if / = 0.1, 10 and 100. Also shown 
(dashed lines) are estimates of the amount of heating that can be 
provided within the cooling radius by the rmal conduction, a ssum- 
ing the radial temperature distribution of lAllen et al.l feOOll) , and 
different values of the suppression factor relative to the maximal 
Spitzer rate (see Section [8] for details). 



8 IMPLICATIONS FOR COOLING FLOWS 

If radio-loud AGN heating is to exactly balance cooling in 
all systems, then the efficiency factor / must increase by 
a factor of 100-1000 between poor groups and rich clus- 
ters. As discussed in Section 17.21 the factor / includes any 
uncertainty or variation in the conversion from radio lumi- 
nosity to jet kinetic energy, and also in the efficiency with 
which the jet kinetic energy is used to heat the cluster within 
the cooling radius; either of these factors may be dependent 
upon environment. In the former case, any boosting of the 
radio luminosities of sources in richer environments due to 
confinement of the radio lobes would work in the opposite 
sense to that required to explain the trend in Figure ITU The 
efficiency of converting the radio source energy into heat, 
however, offers a more feasible environmental dependence of 
/. In the most massive clusters it seems quite reasonable 
that essentially all of the jet kinetic energy is used to heat 
the intracluster medium within the cooling radius of the 
cluster, since the cooling radius is large and observations 
indicate that shock/sound waves are efficiently produced. 
In smaller systems, however, the radio source heating may 
be more inefficient, either because of inefficient production 
of sound/shock waves or because the inflated bubbles rise 
beyond the cooling radius before much of their energy is 
transferred. 

A lower efficiency of converting AGN energy into heat in 
smaller systems would help towards explaining the trend in 
Figure [14] However, given that r adio-mode heatin g is domi- 
nated by low luminosity sources l|Best et al.| [2006). and such 
low luminosity sources tend to have small physical sizes, it 
seems unlikely that it can account for the required factor 



of 100-1000. If it does not, the conclusion has to be that it 
is the balance between AGN heating and radiative cooling 
which changes. In that case, either radio-loud AGN heat- 
ing exceeds that required to balance cooling in groups (but 
matches in rich clusters; / ~ 10), or it falls short of that re- 
quired in rich clusters (but matches that in groups; / ~ 0.1), 
or both (/ ~ 1). 

The possibility that the radio-mode heating in groups 
exceeds that required to balance cooling out to r coo i is in- 
teresting. The luminosity-temperature relation of groups 
is observed to have a steep slope (L <x T 5 ; e.g. 
iHelsdon fc Ponmanll2000l ). much steeper than the slope of 
2 that would be expected from simple scaling laws for viri- 
alised gas, suggesting that the intra- group medium has also 
been heated by a non-gravitational energy source. The ex- 
cess heating required corresponds to about IkeV per par- 
ticle, and has been referred to as the 'entropy floor' (e.g. 
iPonman et al.l B-999). One possibility is that the intra-group 
gas was pre-heated by supernovae from the e arly episodes 
of star formation in the group ga laxies (e.g. iKaiserl Il99ll ; 
iBalogh et aTll 19991 ; iKav et al.ll2003h . However, to obtain suf- 
ficient energy from supernovae requires an extremely high ef- 
ficiency of supernova feedback, and it is also difficult for such 
models to maintai n consistency with the observed metallic - 
ity of that gas (e.g. lWu et al]|200d ; lKravtsov fc Yepesll2000l ). 
An alternative heating sou rce, such as AGN, is there fore 
favoured (see the review bv lMathews fc Brighentill2003l . for 
a full discussion), and the overheating of the intragroup 
medium suggested h ere (for / > 0.1) may provide this. 
ICroston et alj (|2005l ) found that the gas in galaxy groups 
containing radio-loud AGN is typically hotter than that in 
radio-quiet groups, providing a direct indication that radio- 
sour ce heating can be v ery significant; more recently, how- 
ever, IJetha et all (|2007l ) found no temperature difference in 
the very central regions of groups with and without radio 
sources, indicating that any radio source heating occurs pref- 
erentially at larger radii. 

In order for radio-heating to balance cooling in the 
richest clusters, an efficiency factor / > 10 is required. 
If all of the cavity enthalpy is used to heat the intra- 
cluster medium, and all of that heat is transferred within 
the cooling radius, then a value / ~ 4 is expected. If 
there is further heating from weak shocks, or errors in 
the mechanical luminosity determinations (e.g., significant 
heating from sources with extreme mechanical-to-radio lu- 
minosity ratios), then a still higher value could be ob- 
tained, and it may be possible to reach the efficiency re- 
quired. Alternatively, an additional source of heating may 
be present in the most massive systems. The most likely 
process to provide this is thermal conduction, which will 
work to transport energy from the reservoir of hot gas 
outside the cooling radius down into the cooler cluster 
centre (e.g. iNaravan fc Medvedevl l200ll; IVoigt et alj l2002r , 
IZakamska fc Naravanl 120031 ; IVoigt fc Fabian! 12004 and ref- 
erences therein). The heating rate due to thermal conduc- 
tion wa s origi nally calculated for a pure hydrogen gas by 
ISpitzerl (I1962T I, and on energy grounds alone thermal con- 
duction at the maximal Spitzer rate can be sufficient to pro- 
vide all of the heating necessary to balance cool ing in some 
clusters (|Voigt fc Fabian! 12004 IPope et al.| [2006). However, 
in the presence of magnetic fields, thermal conduction will 
be suppressed; the amount of this suppression remains an 
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open question, but iNaravan fc Medvedevl (120011 ) argue that 
in the presence of turbulent magnetic fi elds it may be only 
a factor of a few. Fab ian et al.l |2006) support the argu- 
ment that thermal conduction is relatively efficient in the 
inner regions of clusters, as this is necessary to explain why 
the s hocks seen in th e Pers eus cluster are isothermal. Fur- 
ther, iRevnolds etafl (120051 ) argue that viscosity at about 
25% of the Spitzer value will make the radio bubbles in 
the intracluster medium stable against Rayleigh-Taylor and 
Kelvin-Helmholtz instabilities (although magnetic fields of- 
fer another pos s ibility for providin g this required stability; 
iDe Youndl2003l : iKaiser et al.ll2005h . 

Thermal conduction can t herefore play an impo rtant 
role in heating ga laxy clusters. IVoigt fc Fabianl (|2004l ) and 
iPope et al.l l|2006l ) show, however, that thermal conduction 
cannot provide sufficient heating to balance cooling in all 
clusters, and nor does the radial distribution of the heat- 
ing match that required. In particular, the temperature gra- 
dients in the inner regions of the clusters are too shallow 
to provide sufficient inward he at flux at the clus ter cen- 
tres (see also the discussion of iFabian et alj 120061 . for the 
Perseus cluster). A combined model of AGN heating in the 
inner regions of the cluster, and heating by thermal con- 
duction closer to the cooling radius, would solve this prob- 
lem, howev er. Such a double-heating mod el was first de- 
veloped by iRuszkowski fc Begelmanl ([2002) , and has been 
shown to provide good agreemen t with observed cluster 
properties (see also iBriiggerJ 120031: iHoeft fc Briiggenl [2004 ; 
iRovchowdhurv et al.ll2005l ; iFuiita fc Suzukill2005l ). The rel- 
ative importance of ther mal conduction should al so increase 
in more massive clusters (|Hoeft fc Briiggenl liool ), as would 
be required to account for the decreasing importance of AGN 
hea ting. 

lAllen et alj (|200ll ) studied the temperature profiles of 6 
relaxed clusters observed with Chandra, and found that all 
could be reasonable well-fitted using a universal tempera- 
ture profile. Assuming that this profile is appropriate for all 
cluster^ then the total flow of energy to within the cool- 
ing radius can be approximated, for different assumptions 
of the Spitzer suppression factor, /spitz- These calculations 
are detailed in Appendix [A] and the results are displayed 
on Figure [T4l these confirm that for values of /s P itz ~ 0.25 
thermal conduction could provide a large fraction of the nec- 
essary heating in the most massive clusters, but has little 
importance for clusters with velocity dispersions below 600- 
800 kms" 1 . 

Finally, it is intriguing that for the expected values of 
/ ~ 1 — 4, the combination of thermal conduction (with 
/spitz ~ 0.25) plus AGN heating approximately balances 
cooling within the cooling radius for all systems larger than 
about 400 kms -1 . In this scenario, at lower velocity disper- 
sions the BCG radio-heating would exceed the radiative 
cooling losses within the cooling radius, and for systems 
with velocity dispersions below about 300 kms -1 it would 
exceed the radiative cooling losses of the entire intragroup 
or intracluster medium (not just that within the cooling ra- 



9 The clusters studied by I Allen et al 1 jut ii t ll i had velocity disper- 
sions ranging from 800 to nearly 1500 kms -1 , and so the validity 
of these calculations for systems with velocity dispersions much 
below 800 kms -1 is less certain. 



dius). At this point, the radio source activity would lead to 
heating of the gas, possibly causing some of it to become un- 
bound. This change in the gas properties may then account 
for the change in the luminosity-temperature relation be- 
tween intragroup and intracluster gas that occurs at about 
300-400 kms -1 velocity dispersion. 



9 CONCLUSIONS 

The conclusions of this work can be summarised as follows: 

• Brightest group or cluster galaxies of all stellar masses 
are more likely to host radio-loud AGN than other galax- 
ies of the same mass. The probability of a BCG hosting a 
radio-loud AGN scales roughly linearly with stellar mass, 
which is shallower than the relation for all galaxies. BCGs 
are an order of magnitude more likely to host radio-loud 
AGN than other galaxies at stellar masses below 10 M©, 
but less than a factor of two more likely at stellar masses 
above 5 x 10 n M Q . 

• The distribution of radio luminosities of BCGs is in- 
dependent of the mass of the BCG, and the same as that 
determined for other galaxies. Only the normalisation of the 
radio luminosity function changes, not its shape. 

• The fraction of BCGs of a given stellar mass which host 
radio-loud AGN, and the distribution of their radio lumi- 
nosities, do not depend strongly on the velocity dispersion 
of their surrounding group or cluster. 

• Group and cluster galaxies other than the BCG show 
the same radio properties as those of field galaxies, except 
for those within 0.2r2oo of the centre of the system, which 
show a boosted likelihood of being radio-loud. 

• The fraction of galaxies with emission-line AGN activ- 
ity is smaller in all group and cluster galaxies within r2oo 
than in field galaxies of the same mass. The suppression 
increases with decreasing cluster-centric radius, reaching a 
factor of 2-3 in the centre. With the exception of BCGs at 
the centre of strong cooling flow clusters, BCGs show the 
same emission-line AGN fractions as other galaxies of the 
same mass near the centre of the group / c luster. 

• These results support the argument of (|Best et al.ll2005l ) 
that low-luminosity radio-loud AGN activity and emission- 
line AGN activity are independent physical phenomenon. 
It is argued that the radio-loud activity is associated with 
the cooling of gas from the hot envelopes of elliptical galax- 
ies and, in the case of central cluster galaxies, also from 
the intracluster medium. Accretion of hot gas from a strong 
cooling flow is able to explain both the boosted likelihood 
of BCGs hosting radio-loud AGN, and the different slopes 
of the mass-dependencies of the radio-AGN fractions for 
BCGs and other galaxies. 

• Within the cooling radius or a group / cluster, the me- 
chanical heating output associated with BCG radio-AGN 
activity exceeds that of all other cluster galaxies combined. 

• Either the mechanical-to-radio luminosity ratio or the 
efficiency of converting the mechanical energy of the radio 
source into heating the intracluster medium must be a fac- 
tor 100-1000 higher in rich clusters than in poor groups 
in order that radio-AGN heating balances radiative cool- 
ing for systems of all masses. If not, then radio-loud AGN 
heating either overcompensates the radiative cooling losses 
in galaxy groups, providing an explanation for the entropy 
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floor, and / or falls short of providing enough heat to coun- 
terbalance cooling in the richest clusters. Thermal conduc- 
tion could provide the extra energy required in the richest 
clusters. 
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APPENDIX A: THE HEATING RATE DUE TO 
THERMAL CONDUCTIVITY 

The volume heating rate due to thermal conductivity from 
an approximately infinite heat bath outside a radius r is 
given by 

« 

where k is th e thermal con ductivity of the gas and T is its 
temperature. ISpitzerl (|l962h calculated the thermal conduc- 
tivity for a pure hydrogen gas to be: 

rp5/2 

ks = L84 x iO" 10 ^- Wm^K" 1 (A2) 

where In A is the Coulomb logarithm. In the presence of 
magnetic fields, however, the thermal conductivity will be 
lower. This is usually accounted for by incorporating a factor 
/spitz, such that 

K = Ks/Spitz- (A3) 

The Coulomb logarithm in Equation lA2l depends weakly 
on the gas density and temperature, but to a good approx- 
imation can be treated as a constant with In A w 37 for 
typical conditions in clusters. Hence, 

k = k T 5/2 , (A4) 

whe re Kg fa 5 x 1 0~ 12 ./ "s P itz • 

lAllen et al.l (|200ll ) showed that the Chandra data of 
massive relaxed clusters could be fitted by a 'universal' tem- 
perature profile: 
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T — T2500 



Tb+Ti 



1 + (l/lc)' 



(A5) 



where T = 0.40, T x = 0.61, x c = 0.087 and rj = 1.9 are 
dimensionless constants that they fitted to their data, T2500 
is the temperature of the cluster gas at the radius r2soo 
(within which the mass density is 2500 times the critical 
density of the universe) , and the dimension less parameter x 
is given by x = r/r^snn. [Allen et all l|200ll ) assume r25oo ~ 
0.3r2oo- 

Substituting Equations IA4I and IA5I for k and T into 
Equation IA1I gives: 



T () + Ti 1 



1 + (x/x c 0) v 



To + iTo+Ti) 



x \ r > 



2(1-,) (T +T 1 )(J) 



+ (4T + (2 + 7jj) Ti) (— ) " + 2(1 + 7,) To 



2^ 



(A6) 



The total heating power due to thermal conduction in- 
side a radius r is then roughly 47r?- 3 e(r)/3. Setting r = r- coo i, 



and using r coo i sa 67(<r v /500 kms~ ) kpc from Section [?~3l 
then gives an approximation for the flow of heat to within 
the cooling radius. 

In order to calculate this, r2oo and T2500 are required. 
T2oo can be estimated from the velocity dispersion as r2oo ~ 
1230(cr v /500kms~ 1 )kpc (from iFinn et al.fl2005l . appropri- 
ate for redshift zero and the presently adopted cosmology). 
This then gives 



0.18 



( * V 

1 500 kms" 1 / 



(A7) 



Similarly, T2500 depends on the velocity dispersion, and 
can be estimated from the T — a v relation for clusters 
jXue fc Wull2000h : 



r -° - 2 ' 4x 10? (sooS^r) K - 



(A8) 



Substituting for x and T2500 into Equation IA6I allows 
the total heating rate within the cooling radius due to ther- 
mal conduction to be estimated, for different values of /s P itz- 
Note that strictly these are only valid for systems with ve- 
locity dispersions a v > 800km s -1 , sin ce the 'universal' tem- 
perature profile of I Allen et al.l (|200lh was only determined 
from massive clusters and may not fit lower mass systems. 
Nevertheless, it provides at least an indication of the likely 
importance of thermal conduction within systems of differ- 
ent masses, as shown on Figure [141 
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